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Abstract 

We perform a detailed study of the consistency between different sets of polarized deep inelastic 
scattering data and theory, from the standpoint of a next to leading order QCD global analysis, 
and following the criteria proposed by Collins and Pumplin. In face of recent suggestions that chal- 
lenge the usual assumption about parent parton spin independence of unpolarized fragmentation 
functions, we specially focus on polarized semi-inclusive data. 
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I. INTRODUCTION 



Ever since the renewed burst of interest in the spin structure of the nucleon, triggered 
fifteen years ago by the measurement of the proton spin dependent deep inelastic scattering 
structure function gf by the EMC experiment polarized deep inelastic scattering has 
evolved into a very prolific field for combined theoretical and experimental efforts [2]. As 
result of this activity, an increasingly precise QCD improved partonic description of polar- 
ized nucleons has emerged, phrased in terms of more refined and exhaustive extractions of 
polarized parton distributions. 

Indeed, the number of QCD analyses of polarized inclusive deep inelastic scattering (DIS) 
data at next to leading order (NLO) accuracy has rapidly grown in the few last years ^. 
In some cases these studies also take into account polarized semi-inclusive deep inelastic 
scattering (SIDIS) data 0,0,0]. The usual outcome from these analyses are different sets 
of parton distributions which reproduce fairly well most data sets within the quoted errors, 
with overall values pretty close to the number N of degrees of freedom (d.o.f ) as required 
by the familiar 'hypothesis-testing criterion'. 

Although the apparent fulfillment of this requirement, the resulting parton distributions 
not only suffer from large uncertainties, but in some cases small subsets of inclusive data 
are badly fitted. This facts point to the need of a more stringent criterion for assessing the 
goodness of a particular fit, and also the compatibility between data sets. 

In the case of polarized semi-inclusive data, the motivation for a more careful analysis 
is twofold. For one side, inclusive data do not provide enough information for a complete 
separation of the quark and antiquark distributions of different flavors, so the analyses must 
rely either on SIDIS data or on some external assumption. 

From another side, the standard procedure used to analyze semi-inclusive data, in terms 
of polarized parton densities and unpolarized fragmentation functions, has been criticized 
questioning, in flrst place, the accuracy of the present generation of fragmentation functions 
a. suggesting that the usual assumpt.ous about sp,u indepeudence of fragmentation process 
may not hold |8[, and also addressing the issue of target fragmentation contributions j9||. 

Regarding the flrst issue, recently there has been an increasing amount activity around 
the issue of the polarization of sea quarks in the proton Q]. This interest has been driven in 
part by the conflrmation of the isospin symmetry breaking at sea quark level in unpolarized 
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DIS by the level of refinement attained in QCD global analysis of spin dependent data 
I and by the expected precision of forthcoming polarized semi-inclusive measurements, 
polarized Drell Yan dilepton production, and prompt photon production jl2| . 

Moreover, in spite of the encouraging achievements in polarized semi-inclusive measure- 
ments in the few last years 

HQ, 

which have shown to be in good agreement with inclusive 



data [J| , the available data fail to constrain unambiguously the polarization of sea quarks in 
the proton when included in global QCD fits performed following just the hypothesis testing 
criterion P]. 

In face of this situation, it is worthwhile assessing the goodness of global fits to polarized 
data, involving both inclusive semi-inclusive measurements, with more stringent criteria, in 
the line of what has been proposed in reference and applied to LO unpolarized sets of 
data in The main point of this approach is to apply the parameter fitting criterion to 
subsets of data in the global fit, assessing in this way features as the degree of compatibility 
between subsets of data, the impact of specific data in certain parameters of the fit, and 
the overall consistency of the global fit. The occurrence of such an inconsistency may hint, 
besides an unexpected error source in the experiment, the non validity of assumptions or 
other inputs of the theoretical calculation. 

In the following, after defining our conventions and sketching the computation of the 
relevant observables, we present a global fit of polarized inclusive and semi-inclusive data 
exploring the dependence of the overall Xtot of the fit on the respective xf for each particular 
experiment by means of the Lagrange multiplier method Qj. The analysis is performed 
in two stages. In the first one we restrict the analysis to inclusive data, and to those 
parton densities which can be extracted there. In the second stage we add semi-inclusive 
observables and discriminate between valence and sea quark densities of different flavors, 
but at variance with reference P], where semi-inclusive data was only allowed to flx the sea 
quark polarization, here we leave all the distributions free. 

In order to circumvent the extremely time-consuming convolutions integrals characteristic 
of the semi-inclusive observables at NLO, we apply the double Mellin transform approach, 
developed recently in ^. This new approach has produced results in complete agreement 
with those obtained with the convolution method, but obtained considerably faster. 

As result of our analysis, we flnd that the parameter fltting criterion shows a reasonably 
good level of internal consistency between the inclusive measurements. The addition of SIDIS 
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data sets leads to global fits with /d.o.f. values fairly close to unity, however the increase 
of the statistical weight of this data in the fits shows rather different consequences when 
either positively (cr'^+) or negatively {cr^^) charged final state hadron data is considered. 
For data, the increase in weight does not modify significantly those parton distributions 
coming from inclusive sets, suggesting consistency between the different data sets and also 
between data and the theoretical framework. However, the analysis of data shows a 
sizable degree of conflict with the rest of the fit. 



II. FRAMEWORK. 



Throughout the present analysis, we follow the same conventions and definitions for 
the polarized inclusive asymmetries and parton densities as in reference j^. In the totally 
inclusive case, the spin dependent asymmetries are given by jj] 



(1) 



where the inclusive spin-dependent nucleon structure function gi{x, Q"^) can be written at 
NLO as a convolution between polarized parton densities for quarks and gluons, Aqi{x, Q"^) 
and /S.g^x.Q'^), respectively, and coefficient functions ACj(x)|l 

as{Q_ ^ _ 



A9(i,Q') + 



111 z \ z 



,x 



+ AC,(z)A^?(-,Q^ 

z 



(2) 



A more detailed discussion about these coefficient functions and their factorization scheme 
dependence can be found in Ref. 1^. F^{x, Q"^) is the unpolarized nucleon structure function 
that can be written in terms of F^{x, Q^) and R, the ratio of the longitudinal to transverse 
cross section 

Analogously, for the semi-inclusive asymmetries we have: 

J,dzg^^ix,z,Q') 



Ar{x,Q'')\z'. 



(3) 



dzFC\x,z,Q^y 

hi where the superscript h denotes the hadron detected in the final state, and the variable 
z is given by the ratio between the hadron energy and that of the spectators in the target. 



The region over which z is integrated, is determined by kinematical cuts apphed when 
measuring the asymmetries. 

For the spin dependent structure function 5f(^(x, Q^), we use the NLO expression 



q,q L Jx J z y 



X 



X 



g, Q^) + A, (|, g^) AC« (X, Q^)Df (|, 



, (4) 

and in order to avoid the convolution integrals in x and z we switch to moment space in both 
variables as suggested in Q]. In moment space, the convolution integrals reduce products of 
the Mellin moments of the parton densities 

A/f(Q2)^ ['dxx'^-'AMx,Q'), (5) 
Jo 

fragmentation functions 

ADf "(g2) = /" dxx'^-^AD'l{z,Q^) , (6) 
Jo 

and the double Mellin transform of the coefficient functions ACj- ^^(x, 2;, Q^), defined by 



ACg^'""(g2)= f dxx""-^ f dzz"'-'AC^f{x,z,Q' 
Jo Jo 



(7) 



These coefficients can be written as l2]J 



1 

— + 



1 il + m + nY-l 



(m + 1)2 m(m + l)n(n + 1) 

1 1 
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(9) 
(10) 



where we have set the factorization and renormahzation scales to Q^. As usual, Cp — 4/3, 
Tr = 1/2, and 

^.H-E^. (11) 

For u and d quarks plus antiquarks densities at the initial scale Ql — 0.5 GeV^ we propose 

(/\ A-) - N x''''{l- x)^''{l+-fgX^'') 

'B{a, + l,P, + l)+^,B{a, + 6g + l,(3, + iy 

where B{a,(3) is the standard beta function, while for strange quarks plus antiquarks we 



use 



x{As + As) = 2Ns 



B{as + l,Ps + iy 



with a similar parametric form for gluons 



xAg = Ng 



x"s(l -a;)^9 



(12) 



(13) 



'B(ag + l,(3g + iy 

The first moments of the quark densities 6q (Ng) are often related to the hyperon beta 
decay constants F and D through the SU(3) symmetry relations 



Su + 6u-Sd-Sd = Nu- Na 

= F + D = 1.2573 



(14) 



5u + 5u + 5d + 5d-2{5s + 5s) = Nu + Nd-4.Ns 
^3F-D = 0.579. 



(15) 



Under such an assumption, the previous equations would strongly constrain the normaliza- 
tion of the quark densities. However, as we are not interested in forcing flavor symmetry. 
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we leave aside that strong assumption and relax the symmetry relations introducing two pa- 
rameters, esj and esu{3) respectively. These parameters account quantitatively for eventual 
departures from flavor symmetry considerations, including also some uncertainties on the 
low-x behavior, and higher order corrections, 

N^-N, = {F + D)il + eBj) (16) 



N^ + Nd- 4iV, = (3F -D)il + e5c/(3)), (17) 

and we take them as a measure of the degree of fulfillment of the Bjorken sum rule 2^ and 
the SU (3) symmetry. 

Equations (fTB|) and ()17|) allow to write the normalization of the three quark flavors in 
terms of A^^, ^bj, and esu{3)- Notice that no constraints have been imposed on the breaking 
parameters since we expect them to be fixed by data. The remaining parameters are con- 
strained in such a way that positivity with respect to GRV98 parton distributions is fulfilled. 
These last parameterizations are used in order to compute the denominators of equations 
(1) and (3). This is particularly relevant at large x, and since no polarized data is available 
in that kinematical region, we directly fix the parameters = 3.2, (3^ = 4.05 and (3g = 6 for 
the NLO sets in agreement with GRV98. Consistently with the choice for the unpolarized 
parton distributions, we use the values of Aqcd given in Ref.|22l] to compute ag at NLO. 

As antiquark densities we take 

for u and d quarks, and we assume s = s since the possibility of discrimination in the s 
sector is beyond the precision of the data (as in the unpolarized case). 

Fragmentation functions are taken from [24! and we also use the fiavor separation criterion 



proposed there, which have shown to be in agreement with most recent analysis 

The data sets analyzed include only points with > 1 GeV^ listed in Table 1, and total- 
ing 137, 118, and 34 points, from proton, deuteron, and helium targets respectively, in the 
inclusive stage, plus 42, 24, and 18, from proton, deuteron, and helium targets respectively, 
in the second stage. 
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TABLE I: Inclusive and Semi- inclusive Data used in the fit. 



Collaboration 


Target 


Final state 


# points 


Refs. 


EMC 


proton 


inclusive 


10 


yj 


SMC 


proton, deuteron 


inclusive 


12, 12 


[27] 


E-143 


proton, deuteron 


inclusive 


82, 82 


[28] 


E-155 


proton, deuteron 


inclusive 


24, 24 


[29] 


Hermes 


proton, helium 


inclusive 


9, 9 


M 


E-142 


helium 


inclusive 


8 


m 


F,-1 ^4 


nfil m TTi 
iidi Lliii 


1 n r^l 1 1 Q.^'\T(^ 

llLKjl U-Ol V t3 


1 7 




SMC 


proton, deuteron 


h+, h- 


24, 24 


[13] 


Hermes 


proton, helium 


h+, h- 


18, 18 


[14] 


Regarding the fittin 


g procedure, this is 


done minimizing a 


function /(Ai, A2, . 


., A„) defined 



as 

/(Ai,A2,...,A„)=^A,x^ (19) 

i 

where the sum runs over the data sets or experiments i included in the fit. The parameters 
Aj are the Lagrange multiphers defined for each data set and which can be varied in order 
to produce different fits where the relative weight of a given set of data has been modified. 
Xi is the contribution to the overall coming from the data set i. 

As it is well known, there are various alternatives for calculating these last contributions 
3^ . The most simple and commonly used in fits to polarized data is adding the reported 
statistical and systematic errors in quadrature. This ignores the correlations between data 
points from the same measurements, but in many cases the full correlation matrices are not 
available. 

As inclusive and semi-inclusive data are strongly correlated, and the correlation matrices 
are available, we have taken into account them, analyzing only the inclusive data for SMC 
and Hermes corresponding to 'averaged' bins, for which the correlation matrices are defined. 
In the semi-inclusive case we only consider in the fit the most precise data concerning 
the production of charged ± hadrons (without identifying pions, kaons, or other particles 
individually) . 
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TABLE II: xlW and xjndusiveW different data sets. 





EMC-p 


4.49 


SMC-p 


3.84 


SMC-d 


14.44 


E-143-P 


60.70 


E- 143-d 


83.38 


E-142-n 


4.68 


E-155-P 


17.15 


E-155-d 


17.10 


E-154-n 


6.91 


Hermes-p 


5.15 


Hermes-He 


6.76 



A = 1 A = 20 A = 100 

Xlnc 



224.59 





X/nc 


xf 


X/nc 


4.17 


227.85 


3.46 


262.71 


3.32 


225.21 


3.20 


226.71 


12.81 


245.92 


12.10 


264.77 


60.22 


224.83 


60.21 


224.95 


80.84 


235.32 


80.17 


248.08 


2.40 


233.52 


1.51 


251.98 


16.24 


227.03 


16.23 


233.55 


16.89 


225.65 


14.98 


277.42 


4.19 


229.16 


3.97 


232.02 


4.30 


225.60 


2.83 


234.42 


6.39 


228.31 


5.67 


248.09 



Regarding the normalization uncertainties, in reference j5[ it was found that Hermes data 
was systematically below the best global fits, but allowing a floating normalization factor 
to it in order to account for the relative normalization uncertainties, the values were 
considerably improved. In the following we set allow this factor to be fixed by the global fit, 
finding the best fits for 12% correction. 



III. INCLUSIVE DATA 



In this section we present results from the first stage of our analysis, which only deals with 
inclusive data. It is customary in NLO fits to inclusive data to present several sets of parton 
distributions with different alternatives for the features that are poorly constrained by the 
data, such as the gluon or strange sea quark polarization. As we are mainly concerned 
in assessing the goodness of the fit between theory and data and the degree of internal 
consistency of the data, rather than covering the different scenarios for parton densities, we 
just explore the more favored scenario of reference J], which was labeled as 'set i' and has 
moderate gluon polarization. 
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FIG. 1: xfW Against xjndusiveW f^''" inclusive data sets. 

In Figure H we show the outcome of different global fits to all the available inclusive data. 
The plot give the variation of the xl value of each experiment against the total value of 
the fit ixjuciusive) ■ '^^^ fi^s^ point to the left of each curve (A, = 1) correspond to standard 
fits where no extra weight was given to none of the data sets. The parameters for this fit 
are presented in Table Al in the appendix. Along the curves, the subsequent points come 
from fits where, following the Lagrange multiplier method explained in the previous section, 
increasing values of Aj have been given to a specific set of data, while keeping the others 
parameters \j equal to 1. 

The normal expectation in a good fit to data sets that individually only determine a 
small fraction of the parameters, is a monotonic decrease of a few units in the xl of 
subset which has been subject to the increase in weight while xlndusive varies in the range 

± \/2N. One would also expect Xi fo approach a saturation point within an increase of 
xjnciusive Smaller than V2N. 

As it can be seen in Figure ^ in all the cases, the curves show the features expected for 
consistent subsets of data, each one able to fix a limited number of parameters. The xlW 
values of each experiment i are shown in Table 2. 

Notice that the initial fall in the xlW values with a very mild variation in xjndusive iHus- 
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FIG. 2: Fit to inclusive data (X = 1) together with the Hermes-p and E155-p driven fits with 
A = 50. 

trates the situation where a subset of data would be able to effectively fix some parameters 
in the fit, in general agreement with the remaining data, if it statistical significance were 
increased. This is clearly the situation of neutron and deuteron target data. 

The standard Aj = 1 fit (solid line) can be seen in Figure |21 together with the Hermes- 
p driven fit for Xuermes = 100 (dashes) and that for E155 with A£;i55 = 100 (dots) for 
comparison. As it can be noticed, the changes in the asymmetries due to the extra weight 
in these subsets of data are almost negligible, as it can be expected from the moderate 
variations they produce in XindusiveW- For these reasons we conclude that inclusive data 
is internally consistent and in excellent agreement with theory. 

IV. SEMI-INCLUSIVE DATA 

In this section we focus on the consequences of including SIDIS data in the global fits. 
These data allow in principle to discriminate between light sea quark flavors, so the corre- 
sponding parton distributions are now parameterized and fitted. As we have anticipated, the 
inclusion of SIDIS data leads to acceptable fits according to the hypothesis testing criterion 
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FIG. 3: XiW against XtotW Z*^*" semi-inclusive data sets 

, however the small statistical impact of SIDIS data relative to DIS data hinders definite 
conclusions on sea quark distributions. 

A closer examination allowing for example increased weights in the different SIDIS data 
subsets, as we did for inclusive sets reveals some interesting features, as shown in Figure 3. 
Again, we have applied a 12 % normalization factor to Hermes data, and in order to simplify 
the analysis, we consider SMC and Hermes proton target data together. Identical results 
are obtained if the data is discriminated for each experiment. 

Although for most SIDIS data sets, xli^) reach their respective saturation values within 
a \/2N shift in Xtot required for overall consistency, data coming from positively charged 
hadropro duct ion on proton targets [cTp'^) seems to lay in the borderline, with a fall of several 
units and a saturation value outside de ^/2N range. SIDIS data is mainly dominated by 
proton target data but while a^^ driven fits lead to asymmetries in agreement with other 
data sets, those driven by a^^ show an increasing disagreement with them, as can be seen 
in Figure m 

The inconsistency between these two data sets can also be seen in the parton distributions 
they produce, as shown in Figure El (a) and (b) , for 0"^" and 0"^"^ driven fits at = 5.0 
GeV^, respectively. Notice that (An + Au) and (Ad + Ad) distributions, that should be fixed 
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10"^ 10"' 10"' 1 



X X 

FIG. 4: Fit to SIDIS data (Xi = 1) together with the 0"^+ driven fits. 

mainly by inclusive data, have minor, although not negligible, changes in the fits driven by 
either data sets, showing the degree of consistency between inclusive and semi-inclusive data. 
However, sea quark distributions depend strongly on which data set has received additional 
weight. The discrepancy is particularly strong for Am, which even change sign. 

An interesting feature of SIDIS data is that for all subsets, the SIDIS-driven fits exceeds 
the \/2N range for Xtot with values of A considerably smaller than the ones typical of inclusive 
data, as Figure El shows. While inclusive data sets allows Aj values of 100 or more with a 
few units change in Xtot^ SIDIS data exceeds the allowed range for A > 20 or even less in 
the case of a^~^. In other words, at variance with what happens to DIS data, fits forced to 
reproduce SIDIS data lead to considerably poor global fits. This can be interpreted as a 
weaker level of consistency in the analysis of SIDIS data than in the inclusive case. In Table 
3 we present xf values obtained for each data set in the standard fit (Aj = 1) and increasing 
the weights. 

In addition to the issue of the accuracy of the data, the analysis of SIDIS experiments relay 
also in our knowledge of unpolarized fragmentation functions. Although most of the uncer- 
tainties coming from these functions cancel out when computing asymmetries, the weaker 
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FIG. 5: Parton densities from (a) a'p and (6) £7^+ driven fits at = 5.0 GeV^. 
TABLE III: and xL(^) for different SIDIS data sets. 







A = 1 




A = 10 




A = 20 




xl 


Xtot 




Xtot 




Xtot 




25.61 




21.30 


360.96 


20.94 


450.20 




26.93 




25.57 


332.89 


25.53 


340.10 




5.18 


319.69 


3.99 


385.29 


3.72 


504.24 




14.71 




12.80 


328.27 


12.73 


373.58 


^he 


5.39 




5.36 


319.79 


5.25 


329.17 




6.69 




6.46 


338.18 


5.71 


409.15 



degree of consistency, and odd behavior of cr regardless one uses proton or deuterium 
targets, and in different experiments, may hint a failure in the extraction of fragmentation 
functions, particularly in the more troublesome discrimination between those for positive 
and negative final state hadrons. 
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V. CONCLUSIONS 



We have assessed the internal consistency in a NLO analysis of different sets of polarized 
DIS and SIDIS data. For the inclusive data, the agreement shown between theory and data 
from the stand point of the hypothesis testing criterion, is confirmed when a detailed analysis 
using a variant of the parameter fitting criterion is performed. For SIDIS data the level of 
consistency is considerably weaker, particularly in the case of positively charged final state 
hadron data. Even though this kind of analysis can not establish whereas either data or 
some particular ingredient in the theoretical approach is responsible for the discrepancy, the 
unexpected features found in the analysis of a^'^ data, for different targets and coming from 
different experiments, hint an inaccuracy in the separation between positively and negatively 
charged hadron fragmentation functions. 
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VI. APPENDIX: PARAMETERS OF THE FIT 



We present here the parameters of the fit with A = 1, for the first stage where we only 
deal with the inclusive data {DIS Fit) and the second stage, where the semi-inclusive data 
were included {SI DIS Fit). 



Table A I: Parameters for the lambda = 1 fit. 



Parameter 


DIS Fit 


SI DIS Fit 




-0.002 


-0.004 




0.085 


O.OScS 


OLu 


0.858 


0.858 


fin 


3.200 


3.200 


lu 


14.929 


14.969 




1.004 


1.005 


OLd 


0.434 


0.433 


Pd 


4.050 


4.050 


Id 


13.888 


13.939 


Sd 


1.651 


1.651 


Ns 


-0.074 


-0.075 


as 


2.500 


2.491 


(3s 


10.000 


10.000 




0.239 


0.238 


ag 


1.499 


1.499 


Pa 


6.000 


6.000 


Nu 




-0.014 






2.311 






7.646 






0.014 






2.315 


Pd 




7.646 
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